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Abstract 

The focus of this study was to reconstruct a phylogenetic hypothesis for the moth subfamily Arctiinae (tiger moths, woolly 
bears) to Investigate the evolution of larval and adult pharmacophagy of pyrrollzldlne alkaloids (PAs) and the pathway to PA 
chemical specialization in Arctiinae. Pharmacophagy, collection of chemicals for non-nutritive purposes, is well documented 
in many species, including the model species Utetheisa omatrix L. A total of 86 exemplar ingroup species representing tiger 
moth tribes and subtribes (68 genera) and nine outgroup species were selected. Ingroup species included the most species- 
rich generic groups to represent the diversity of host-plant associations and pharmacophagous behaviors found throughout 
Arctiinae. Up to nine genetic markers were sequenced: one mitochondrial (COI barcode region), one nuclear rRNA (D2 
region, 28S rRNA), and seven nuclear protein-coding gene fragments: elongation factor 1-a protein, wingless, ribosomal 
protein subunit S5, carbamoylphosphate synthase domain regions, glyceraldehyde-3-phosphate dehydrogenase, isocitrate 
dehydrogenase and cytosollc malate dehydrogenase. A total of 6984 bp was obtained for most species. These data were 
analyzed using model-based phylogenetic methods: maximum likelihood (ML) and Bayesian inference (Bl). Ancestral 
pharmacophagous behaviors and obligate PA associations were reconstructed using the resulting Bayes topology and 
Reconstructing Ancestral States in Phylogenies (RASP) software. Our results corroborate earlier studies on the evolution of 
adult pharmacophagous behaviors, suggesting that this behavior arose multiple times and Is concentrated in the 
phaegopterine-euchromline-ctenuchlne clade (PEC). Our results suggest that PA specialization may have arisen early In the 
phylogeny of the subfamily and that facultative larval pharmacophagous behaviors are the derived condition. 
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Introduction 

Arctiinae (tiger moths and woolly bears) are a charismatic moth 
lineage with a complex evolutionary relationship with plant and 
fungal chemistries. Comprised of approximately 11,000 species 
[1], this cosmopolitan group is well known among ecologists and 
evolutionary biologists for the evolution of bright coloration and 
spectacular adult mimicry of wasps, beetles and unpalatable moths 
and butterflies (Figure 1). Both aposematic adults (Figure lA) and 
larvae (Figure IB-D) typically harbor endogenous biogenic 
amines, like histamines, which are often supplemented with 
secondary compounds acquired from larval hosts or through 
pharmacophagy - feeding on plants to obtain chemicals rather than 
nutrients [2]. Some larvae become pharmacophagous when 
parasitized, and this "self-medication" improves survivorship [3]. 
Adult pharmacophagy is linked to acquiring chemical constituents 
of courtship pheromone and nuptial gifts that improve male 
mating success [2], [4-7]. Adult pharmocophagy has been 
documented in other Lepidoptera, notably Danainae (Nympha- 
lidae) [2]; however, larval pharmacophagy has only been 
documented in Arctiinae. 



Within the mostly phytophagous tribes Arctiini and Amerilini, 
larvae supplement endogenous defenses by sequestering com- 
pounds from their host plants like pyrrolizidine alkaloids (PAs), 
cardiac glycosides (CGs) or iridoid glycosides (IRs) [4], [8-10]. 
These noxious compounds are of great importance because they 
have been implicated in severe human and animal diseases [2]. 
Additionally, through behavioral assays of vertebrate and inver- 
tebrate predators, these three classes of compounds (PAs, CGs, 
and IRs) have been shown to deter predation [11], and of these 
classes of compounds, the acquisition, sequestration and dissem- 
ination of PAs is best understood [12]. 

Larvae of PA specialist species (e.g., Uletheisa omatrix, Tyria 
jacohaeae (L.) [13], Creatonotos gangis (L.)) feed only on PA- 
containing plants and prefer foliage and seed pods that have the 
highest concentration of these compounds [12], [14]. Larval PAs 
are stored in an inactive form for transfer to the adult stage. In 
some adult males, PAs are transformed into danaidol (and 
chemical analogs) for dispersal in courtship pheromones [7]. 
Within the Old World genus Creatonotos, PAs are morphogenic; 
the size and complexity of the male-pheromone dispersing 
structures (coremata) are dependent upon the PA concentrations 
consumed by larvae [15]. In the rattlebox moth, Utetheisa 
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Figure 1. Tiger motKi adults/larvae, diversity. A. Diversity of adult habitus (photos courtesy: Rebecca Simmons); B. Woolly Bear, Pyrrharctia 
Isabella (photo courtesy of Bill Conner); C. Milkweed tussock, Euchaetes egle (photo courtesy of Rebecca Simmons), D. Rattle box moth larva, Utetheisa 
ornatrix (photo courtesy Nancy Jacobson). 
doi:1 0.1 371/journal.pone.01 01 975.g001 



ornatrix, male mating success depends upon the presence and titer 
of PAs in his pheromones [16-19]. In this species, the males 
transfer PAs to the female in the spermatophore [20], and the 
females gain predator protection for themselves and their eggs 
[21]. In the rattlebox moth, PAs are considered an "honest" good- 
genes signal of the male's ability to acquire and transport PAs [18], 
[22-23]. 

In contrast to this PA-obligate system, where PA acquisition and 
use is linked across the larval and adult stages, several have been 
chemically documented, and many species have host records of 
PA-containing plants. In oligophagous or polyphagous species, 
larvae accept PA-containing plants as part of a broader diet. They 
do not transfer them to the adult stage or actively collect them as 
adults. In some cases, the purpose of PA acquisition appears to be 
medicinal: larvae increase their PA diet when parasitized by 
endoparasitoids [24-25]. Thus, parasitized larvae that self- 
medicate with PAs tend to have greater resistance to parasitoids 
and thus higher survivorship than those that do not [26]. The 
precise mechanism leading to higher survivorship is not under- 
stood. Higher titers of PAs have also been shown to delay larval 
development and incur other fitness tradeoffs [26-27]. This 
behavioral plasticity of facultative pharmacophagy is thought to 
permit polyphagous larvae to respond to a heterogeneous selective 
environment [3], [25]. Although larvae facultatively acquire PAs, 
they do not store them through metamorphosis and transmit these 
to the pupal and adult stages. Thus, in these cases, the PA 
pharmacophagy is restricted to the larval stage. 

In contrast, there are several species whose pharmacophagy of 
PAs only occurs during the adult stage (Figure 2). Many species 
(typically males) seek out PA-plants to feed at wounds or acquire 
PAs by regurgitating saliva, dissolving PA crystals on withered 
plant leaves and re-ingesting the solution [2], [5-6], [28]. Males of 
Cosmosoma myrodora (Dyar) have deciduous scales impregnated 
with PAs that are discharged from a ventral abdominal pouch to 
cover the female prior to copulation. The female is instantiy 
protected from invertebrate predators by her PA "wedding veil" 
that is released by the male during courtship [29-30]. Many of the 



species that exhibit adult pharmacophagy occur in the Neotropics, 
with the exception of the Old World genus Amerila. 

Despite numerous ecological and behavioral studies, our 
understanding of arctiine phylogeny and how these pharmaco- 
phagous feeding behaviors evolved was limited [10]. A culmina- 
tion of a series of molecular and morphological studies led to the 
recent placement of the family Arctiidae as a subfamily Arctiinae 
of Erebidae [31-32] and all taxonomic ranks of older literature is 
translated to this most recent taxonomy (Table 1) [33]. The only 
comprehensive phylogenetic reconstruction for Arctiinae was 
based on adult and immature morphology from 40 exemplar 
species [34]. The results supported the monophyly of the 
subfamily and three lineages, tribe Arctiini and the traditional 
arrangement of a sister group relationship of Lithosiini and 
Syntomini (Figure 3A). Their tree disagreed with Bendib and 




Figure 2. Halysidota tessellaris (Smith) 
courtesy Bill Conner). 

doi:1 0.1 371/journal.pone.01 01 975.g002 
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Minet [35], who proposed that Syntomini + Arctiini was sister to 
Lithosiini (Figure 3B). The Old World genus Amerila was treated 
as part of the Arctiini (see review: [9]) until 2009, when a fourth 
lineage, Amerilini (= Rhodogastrinae of KiriakoflF [36]) was 
resurrected [37]. A recent molecular study [31] suggested that the 
four lineages should be arranged (Lithosiini(Amerilini(Syntomini + 
Arctiini), but only 1 7 arctiinae species were included in the much 
larger study of Erebidae relationships. 

The focus of this study was to reconstruct a phylogenetic 
hypothesis of Arctiinae with a large sample of species spanning the 
diversity of lineages and taxonomy with nine genetic markers. The 
resulting phylogeny is used to investigate the evolution of 
pharmacophagy and the pathway to chemical specialization in 
Arctiinae. 

Materials and Methods 

Taxon Sampling 

Our goal was to sample the most species-rich generic groups to 
represent the diversity of host plant associations and feeding 
behaviors (e.g., pharmacophagy) found throughout Arctiinae. A 
total of 86 species from 68 genera representing tiger moth tribes 
and subtribes were selected based on previous treatments and 
checklists [31-32], [34-35], [37-39]. We sampled from three 
traditionally recognized tribes (formerly subfamilies): Lithosiini, 
Syntomini and Arctiini. We also included Amerilini and the 
recently proposed Spilosomini [38] . Representatives from former 
tribes (now subtribes) included Arctiina, Callimorphina, Perico- 
pina, Euchromiina and Ctenuchina [1]. Representative lyman- 
triines, aganaines, and calpines were used as outgroups (nine taxa). 
Fieldwork related to the present study (e.g., specimen acquisition) 
did not involve endangered or protected species. Field locations 
where permits were required with relevant contact information are 
as follows: USA: Big Cypress National Preserve, FL: collecting and 
destructive samphng permits were obtained through National 
Parks Services (Jennifer Stafford and Steve Schulze); COSTA 
RICA: permits were obtained through the Organization for 
Tropical Studies (Francisco Campos Rivera); RUSSL\N FEDER- 
ATION: Primorsky Krai, permits obtained through the Far 
Eastern Branch of the Russian Academy of Science (Dr. Vladimir 



Kononeko); AUSTRALIA: Lamington National Park, Queens- 
land: permits obtained through the Queensland Government 
Environmental Protection Agency (Jacqui Brock and Ian Bryant). 
Our exemplars were crossed-checked with their COI sequences 
(barcode region) in Barcode of Life Data System (BOLD) [40]. A 
synopsis of the classification, gene sampling and Genbank 
accession numbers for species included in the study is provided 
in Table SI. 

DNA Markers, Extraction, PGR and Sequencing 

DNA extractions were performed on 1—2 legs (or in some cases 
abdominal tissue) using the QIAmp Micro DNA Extraction kit 
(Qiagen, Valencia, CA). The following nine genetic markers were 
amplified based on previous studies demonstrating their usefulness 
in resolving evolutionary relationships among lepidopteran insects 
above and below the family level [41-42]: (1) the D2 region of the 
28S ribosomal subunit sequence, (2) the barcode region of the 
COI gene, (3) elongation factor l-0( protein (EFl-a), (4) ribosomal 
protein subunit S5 (RpS5), (5) wingless (W'GS), (6) carbamoylpho- 
sphate synthase domain protein (CAD) (7), glyceraldehyde-3- 
phosphate dehydrogenase (GAPDH) (8), isocitrate dehydrogenase 
(IDH) and (9) cytosolic malate dehydrogenase (MDH). A total of 
6984 bp was obtained for most species. Polymerase chain reaction 
(PGR) for the majority of the data in this study was performed 
using ACCUZYME Mix (Bioline) using the following modified 
recipe: 12.5 i^L ACCUZYME Mix (Bioline), 0.5 |xL Taq poly- 
merase (Bioline), 1.0 |J.L of each primer (10 |iM), 1.0 |J,L DNA 
template, and 9.0 |xL ddH20. PGR amplification of 28S used 
primer pair 28S-F1 and 28S-R1 [42] and the author's cycling 
profile, except touchdown PGR was omitted. Amplification of 
COI was performed with primers Lep-Fl and Lep-Rl [43]. The 
cycling profile consisted of 3 min at 94°C, 33 cycles of 1 min at 
94°C, 1 min at 48°C, and 1 min at 72°C, and a final extension 
period of 5 min at 72°C. PGR of EFl-ot was achieved using the 
primer pair rcM4 and M46- 1 with the modification to M46- 1 that 
the fifth base from the 3' was made degenerate: G^ (AG) [44]. 
Thermocycling conditions were 5 min at 95°G, 40 cycles of 94°C 
for 30 sec, 50°C for 30 sec, and 72°C for 90 sec, with a final 
extension at 72°G for 10 min. Amplification of RpS5, WGS, 
GAD, GAPDH, IDH and MDH foUowed previously described 



Table 1 . Examples of recent taxonomic treatments of the former lepidoptera 


n family Arctiidae within the superfamily Noctuoidea. 




Bendib & Minet 1998 


Kitching & Rawlins 1998; Jacobson & Weller 2002 


van Nieukerken et al. 201 1 (and references therein) 


Noctuoidea 


Noctuoidea 


Noctuoidea 


Erebidae 




Lymantriidae 


Lymantriinae 


Arctiidae 


Arctiidae 


Arctiinae 


Amerilinae 


Lithosiinae 


Lithosiinae 


Lithosiini 


Syntominae 


Syntominae 


Syntomini 


Arctiinae 


Arctiinae 


Arctiini 




Arctiini 


Arctiina 




Callimorphini 


Callimorphina 




Pericopini 


Pericopina 




Phaegopterini 


Phaegopterina 




Ctenuchini 


Ctenuchina 




Euchromiini 


Euchromiina 



doi:10.1371/journal.pone.0101975.t001 
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Figure 3. Competing Kiypotheses for subfamily reiationsKiips in tiger moths based on a) Jacobson and Weller 2002 and b) Bendib 
and Minet 1998. 

doi:1 0.1 371 /journal.pone.01 01 975.g003 



protocols [31]. To facilitate easier high-throughput sequencing, all 
primer pairs were appended on the 5' end with the universal T7/ 
T3 primer tails. Successfully amplified products were purified 
prior to sequencing either the QIAquick PGR Purification kit 
(QIAGEN) or USB ExoSAP-IT PGR Product Cleanup (AflFyme- 
trix) using the manufacturers' protocols. Purified PGR products 
were sequenced on an ABI 3730x1 (Applied Biosystems) using ABI 
BigDye v. 3.1 Terminator chemistry (PE Applied Biosystems) at 
the BioMedical Genomics Center (BMGG) at either the University 
of Minnesota, St. Paul campus or by the company Macrogen Inc 
(Amsterdam, the Netherlands). 

DNA Sequence Alignment and Phylogenetic Analysis 

All protein-coding sequences were aligned using BioEdit v7.2.5 
[45] and 28S sequences were aligned using MAFFT v.7.017 [46] 
implemented in Geneious R6 v.6.1.6 [47]. DNA sequence data 
was managed and datasets were generated using the web-based 
software VoSeq [48]. Neighbor-joining and Bayesian analyses of 
single-gene alignments were performed to test for rogue sequences. 
If in a single gene analysis, a species placed in a difiFerent 
relationship compared to its placement in the combined, the 
original sequence data were examined to ensure that contamina- 
tion had not occurred. In some instances, the sequence amplified 
but sequenced poorly. If attempts to procure better sequence were 
unsuccessful, questionable data were omitted. 

Maximum likelihood (ML) and Bayesian Inference (BI) analyses 
were carried out on the combined data set of all markers. The 
Bayesian Information Criterion (BIG) using Partition Finder v. 
1.1.1 [49] was apphed to determine the best partitioning scheme 
for codon positions in protein-coding sequences and correspond- 
ing best evolutionary models for the dataset (Table 2). A maximum 
likelihood (ML) analysis was performed using RAxML-HPG 
BlackBox v. 7.4.4 [50] on the GIPRES Web Portal at http:// 
www.phylo.org/portal2/ [51]. Bootstrapping runs (1000 repli- 
cates) to calculate ML nodal support were also performed on the 
GIPRES site using the RAxML-HPG BlackBox utility. A 



partitioned Bayesian analysis of 2.0x10^ generations, with 25% 
of trees discarded as burn-in, was carried out using the MrBayes 
v.3.2.2 [52-53] on XSEDE utility on the GIPRES Web Portal at 
http://www.phylo.org/portal2/ [51]. 

Evolution of PA Acquisition Strategies 

An ancestral state reconstruction analysis of PA acquisition 
strategies was undertaken to elucidate possible feeding shifts and 
origins of pharmacophagy. Known PA acquisition strategies for 
tiger moths were divided into the foUowing functional categories: 
(A) PA-feeding as adult, (B) PA-generalist feeding as larvae, and (G) 
PA-specialist feeding as larvae. PA acquisition strategies (A-G) 
were coded for terminal taxa based on documented cases in the 
primary literature [9]. These categories reflect what is currently 
known about the biochemical mechanisms of pyrrolizidine 
alkaloid sequestration (or lack thereol) in tiger moths that are 
known to be associated with PA plants [12]. Published feeding 
records suggest PA acquisition strategies in tiger moths are 
generally conserved at the genus level [9]. Therefore, in cases 
where published records for a given species were absent, available 
records for congeners were used (approach as in [54]). In some 
cases, species exhibited more than one acquisition strategy and 
thus were assigned more than one state (e.g., Nyctemera and 
Haploa). Although other ingroup taxa likely do acquire (and 
possibly sequester) PA's, detailed documentation of occurrences is 
limited. For convenience, when assessment of PA use was 
unavailable, and terminal species for which genus-level data were 
not available, these species were treated as unknown (NuU). 
Character codings for species in the analysis are included in 
supplemental materials (Table S2). Tables summarizing PA source 
records can be found in Conner and Jordan [9] as well as earlier 
studies [7-8]. 

This data matrix was used to estimate ancestral reconstructions 
for known PA acquisition strategies following previously published 
methods [55] . Analyses were based on the resulting trees from the 
ML analysis of molecular sequence data. We performed a 
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Table 2. Optimal partitioning scheme 


selected by PartltionFlnder vl.0.1 (Lanfear ef al. 2012) using the BIC selection criterion. 




Pdrtition 


Bfist Psrtitioninc) Schems 


Included Nucleotides 


1 


1 ViVI+l+U 


position 1 of CAD, position 1 of MDH 


2 


GTR+G 


position 2 of CAD, position 3 of IDH, position 2 of wingless 


3 


K80+I+G 


position 3 of CAD, position 3 of MDH, position 3 of RpS5, position 3 of wingless 


4 


GTR+I+G 


position 1 of COI, position 3 of Efl a, position 3 and 1 of GAPDH, position 1 and 2 of IDH, position 1 and 2 of RpS5, 
position 1 of wingless, 28S 


5 


K81 uf+l+G 


position 2 of COI 


6 


TrN+l+G 


position 3 of COI, position 2 of Efl a, position 2 of GAPDH, position 2 of MDH 


7 


SYM+I+G 


position 1 of Efl ot 


The model of evolution for each partition, the number of subsets, parameters, and the log-likelihood for the scheme are included. 
doi:10.1371/journal.pone.0101975.t002 



related genera ( = subtribe Callimorphma). Another clade 
comprising Nyctemera + Secusio, was also recovered (BS = 100, 
PP =1). Subtribe Arctiina was recovered with liigh support (BS = 
100, PP = 1). The BI or ML analyses did not recover Spilosomina 
exclusive of Arctiina [38]. 

Witliin Arctiini, a second large clade included the remaining 
subtribes: Pericopina, Phaegopterina, Ctenuchina and Euchro- 
miina, and this clade was weakly supported in the ML while 
strongly supported ui the BI (BS = 19, PP = .98). Within this clade 
(Figure 4), there was strong support for Pericopina (BS =100, PP 
= 1), a small generic cluster oi Melese and Bertholdia (BS = 100, 
PP = 1), and a Eucereon-Ctenucha clade (BS = 98, PP = .99). 
Euchromiina was recovered with moderate support in the ML 
analysis (BS = 75, PP — .87). Phaegopterina was not recovered as 
monophyletic in either analysis, but this result was not surprising 
given prior studies also lacked support for the monophyly of the 
subtribe [34], [59]. Remaining structure was either not well 
supported or conflicted between die two analyses (although nodal 
support was typically weak). 

Ancestral State Reconstructions of PA Acquisition 
Strategies 

The resulting ML topology with the most likely reconstructions 
for PA strategy is shown in Figure 5. Our analysis supports the 
hypothesis of WeUer et al. [10] that the incorporation of PA's into 
the larval or adult diets occurred after the ancestor of lichen 
feeding lithosiines and ancestor of the Arctiini + AmerUini clade 
diverged (Figure 5). The analysis does not reconstruct a single 
origin of larval or adult acquisition nor does it find a progression 
from larval acquisition early in the evolution of tiger moths with 
adult feeding arising later [10]. Rather, adult PA feeding arises 
early with a single origin within AmerUini and at least two 
independent origins in the clade containing Pericopina, Phaegop- 
terina, Euchromiina and Ctenuchina (Figure 5). Species in the 
genera Euchromia and Macroneme are also known to acquire PA's 
during the adult stage but a significant reconstruction of that state 
was not recovered for Euchromiina, likely due to unknown 
behaviors in closely related taxa. 

Obligate larval acquisition and sequestration (e.g. Utetheisa sp. 
and Tyria sp.) precedes the evolution of the generahst feeding 
strategy within CaUimorphina (Figure 5). The callimorphine genus 
Haploa represents a special case of generalist PA feeding within 
this subtribe because larvae begin as PA plant specialists and then 
become polyphagous during later instars [9]. Nyctemera species 
(also CaUimorphina) are generalist PA feeders during both larval 
and aduh stages [2], [9], [14], [60], [61]. Generalist PA feeding 



Bayesian binary MCMC (BBM) analysis implemented in RASP 
(Reconstruct Ancestral State in Phylogenies) [56-57] for 
20,000,000 generations with sampling at every 100 generations. 
This method calculates the frequencies of ancestral distributions, 
ranges or other traits at each node and averages them over all 
sampled trees [56-57]. The fixed JC model and nuU distribution 
were used. The maximum number of feeding categories was set at 
three (with combinations of distributions allowable). The RASP 
analysis was rooted with the same root for the phylogenetic 
estimation, a non-PA feeder C. thalictri (Erebidae; Table S2). 

Results 

Phylogenetic Analyses 

In the analyses of the combined dataset, BI and ML approaches 
produced congruent results. The BI and ML analyses recovered a 
well-supported, monophyletic Arctiinae (bootstrap [BS] = 100, 
posterior probability [PP] = 1) (Figure 4). A clade comprising 
aganaine genera was sister to it with strong support for this 
relationship in both analyses (BS = 100, PP = 1). Within 
Arctiinae, four clades corresponding to the tribes Lithosiini, 
Amerilini, Syntomini and Arctiini were recovered. Tribe Lithosiini 
placed as sister to remaining arctiines (Figure 4) and its monophyly 
was strongly supported (BS =96, PP =1). Within Lithosiini, 
highly structured, well-supported clades were recovered. Genera 
Lyclene and Miltochrista were not recovered as monophyletic 
(Figure 4), but this result was not unexpected given a much larger 
study of Lithosiini that included multiple species from these genera 
[58]. Tribe Amerilini (BS = 100, PP = 1) was recovered and 
placed as sister to the remaining arctiines (BS = 100, PP = 1). 
Support for the clade comprising Arctiini and Syntomini in the BI 
analysis was high (PP = .98) and moderate in the resulting ML tree 
(BS = 76) (Figure 4). However, support for tribe Syntomini was 
strong in both analyses (BS = 100, PP = 1). 

The tribe Arctiini received high support [BS =93, PP = 1], 
however the genus Utetheisa, which has been associated with 
Callimorpha and related genera, differed substantively in its 
placement in the two analyses. In the ML analysis, Utetheisa 
placed as sister to the rest of Arctiini (Figure 4). In the BI analysis, 
Utetheisa placed within the clade comprising pericopines, 
phaegopterines, euchromiine and ctenuchine species (PP = .98) 
with the novel association of genera Pagara and Castrica (PP 
= .90) (Figure SI). This unexpected result is discussed below. 

With the exception of Utetheisa, the two analyses agree on the 
remaining large clades and many of the generic associations. 
CaUimorphina grade into Arctiina in both analyses (Figure 4). 
However, a large clade is recovered comprising Callimorpha and 
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Arctornis sp. 

Orgyia aiitiqua 

Lyman Iria monacha 

Nygmki plana 



Lycleiie reticulata 



Barsine sp. 




O 



'Halysidota tessellaris 
'Pseiidohemihyalea labecula 
Pseudohemihyalea edwardsii 
Elysius rabusculum 

Leucanopsis setosa 
Ciirales arizonensis 
Pliaegoplera decrepidoides 
Praemastus albipimcta 
'Ainastus episcotosia 
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Figure 4. Phylogenetic hypothesis for the subfamily Arctiinae (Noctuoidea, Erebidae) based on maximum lilcelihood (IVIL) analysis, 
along with outgroups. Clades representing tribes are colored. Support values (ML Bootstrap/PP posterior probability) are shown next to the 
branches. 

doi:l 0.1 371/journal.pone.01 01 975.g004 



with supplemental pharmacophagy during the immature stages 
also dominates within Arctiina and may have three origins within 
the Arctiina; however, the sampling of species (a fraction of the 
diversity) and knowledge of these behaviors is such that the pattern 
should not be over-interpreted. 

Discussion 

Phylogeny 

The basal position of Lithosiini is expected based on prior 
studies [31]. Interestingly, the enigmatic Amerila (AmerUini) is 
recovered as sister to Arctiini and Syntomini. This genus has been 
previously placed in Phaegopterina [35], [61-62]. Its placement as 
sister to the non-lithosiine tiger moths was first suggested by Zahiri 
et al. [31] in the context of their larger study of Erebidae. They 
had very limited sampling of one Amerila species, three 
syntomines and nine arctiines, but sampled one mitochondrial 
marker (COI) and seven nuclear genes for their erebid study of 
237 species [31]. In their tree, the node that placed Amerilini as 
sister to Syntomini and Arctiini had 52% ML bootstrap support 
compared to 95% support in our results; they did not perform a 
Bayesian analysis. This placement supports the interpretation that 
the hearing organ morphology, a notodontid-type tympanum, is 
independentiy expressed in Amerila, Syntomini, Creatonotos 
(Arctiini, Arctiina) and several phaegopterine genera in Arctiini. 
Finally, the sister relationship of Syntomini and Arctiini supports 
the hypothesis of Bendib and Minet [35] . 

Within Arctiini, the relationships of genera formerly placed in 
CaUimorphina [Callimorpha, Euplagia, Haploa, Nyctemera, 
Secusio, Tyria, and Utetheisa) need comment. This enigmatic 
and relatively small group of about 100 "butterfly-like" species 
placed in eleven genera has a cosmopolitan distribution [1], [62]. 
The brightly colored, diurnal members are popular with collectors 
and many forms, subspecies and subgenera have been described 
[1]. At least two tribes (now subtribes) have been recognized 
including Nyctemerina, comprising two genera Nyctemera and 
Secusio (Galtara) [37] and CaUimorphina, whose membership has 
been highly variable, but centered on Callimorpha and Euplagia 
[62]. That generic pairs, like Nyctemera and Secusio, are supported 
strongly in our analysis and by prior authors [37], [62-63] is not 
surprising, however, the taxonomic usefulness of a subtribe that 
only includes two of eleven possible genera and about 100 species 
is debatable. Regardless, the remaining CaUimorphina is recov- 
ered as a grade of genera subtending Arctiina whose precise 
generic relationships will require additional species sampling to 
illuminate. The instability of the placement of Utetheisa is likely 
due to the effect of missing species diversity in our sample from the 
Austral-Asian and African regions, as well as poor coverage of 
gene regions (only RpS5, wingless, COI and 28S were sequenced 
for this genus, see Table 81). We include one New World Utetheisa 
{U. ornatrix L.) and Old World tropics U. lotrix (Cramer). In an 
unpublished dissertation [63], a morphological study of adults and 
maximum parsimony analysis of the fourteen species found that U. 
ornatrix originated early in the evolution of the genus and was not 
closely related to the Galapagos Ulelheisa sp(;ci(-s. Rather, the 
species U. lotrix and then U. pulchella L. placed in the clade as 
sister to the Galapagos species. Additional Utetheisa species and 
the sister genus Pitasila should be sampled. Taking taxonomic 



action to redefine CaUimorphina is not recommended until these 
relationships are tested with more species. 

The subtribe Arctiina (formerly Arctiini [34], [59]) is strongly 
supported by our data. Although our taxon sampling was not 
designed to specifically test generic limits, our results do reflect 
some of the arctiine relationships proposed by Ferguson [64] . For 
example, Apantesis and Palearctia form a sister grouping 
(Neoarctia-Grammia group), although the relationship did not 
receive strong support in either analysis. The ^'Holomelina" group 
of Ferguson [64], now Virbia (sensu Zaspel and WeUer [65]), was 
strongly supported as monophyletic in both analyses. This clade of 
four Virbia species robustiy placed within the assemblage of genera 
commonly referred to as the Spilosoma Group of Ferguson [64] or 
SpUosomina (formerly SpUosomini; [37]). The taxonomic useftil- 
ness of this subtribe is questionable given the tendency of workers 
in this group of moths is to create higher taxonomic categories to 
accommodate the diverse color patterns. 

A clade of genera placed in Pericopina, Phaegopterina, 
Euchromiina and Ctenuchina is recovered (Figure 4) with no 
ML support (n.s., 19) and with high BI support (PP =.98). The 
association of these genera have been supported by numerous 
morphological studies and faunal works (review [34]), and the low 
ML support is surprising. This clade represents over half the 
species diversity in the subfamUy. We interpret the low ML 
support as reflecting missing lineage diversity. Within this clade, 
the generic groupings reflect general consensus based upon prior 
studies and taxonomic treatments (e.g., [38], [66-68]), although 
none of these nodes are significantiy supported by the ML analysis. 
Our analyses also recover the Halysidota Group of Forbes [69] 
(ML = 73, PP = .90) and a possible "jamming" clade (ML = 100, 
PP = 1) comprised of species [Bertholdia trigona + Melese spp.) 
known to disrupt echolocation by their bat predators [70-7 1]. The 
analyses recover the Euchaetes Group of Forbes [69] with the 
inclusion of Pagara in the ML analysis (BS = 26) and its exclusion 
in the BI analysis (PP = .91). The placement oi Pagara as sister to 
Castrica (Figure SI) in the BI analysis is a novel association. Given 
the low ML support, we do not place confidence in this alternate 
arrangement. 

Evolution of PA Acquisition Strategies 

Our results suggest numerous shifts to PA pharmacophagy as 
larvae, as adults and once for both life stages within Arctiini. The 
first occurrence of PA adtUt pharmacophagy is within the tribe 
AmerUini (Figure 5). Many species of the Old World genus 
Amerila are known to be pharmacophagous as adults [72], but 
their larval feeding behaviors are unknown. Syntomini and 
Dysschema are not associated with PA feeding as either larvae or 
adtUts. 

Within the grade of CaUimorphina genera, the evolution of 
obUgatory PA larval specialists occurs in Utetheisa, whose larvae 
specialize on PA plants in the family Fabaceae (e.g. Crotalaria) and 
in Tyria jacoheae, a known biological control agent whose larvae 
specialize on Asteraceae (Senecio spp.). The clade comprising 
Callimorpha + {Euplagia + Haploa) shows a pattern of larval 
pharmacophagy (Figure 5) with Haploa being an exception. This 
genus contains some species whose larvae are PA specialists during 
the early instar stages and then shift to a generalist, pharmaco- 
phagous strategy during later instar stages [73], or after 
hibernation [13]. It is unclear whether this pattern of PA 
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Figure 5. Summary of Bayesian ancestral state reconstruction 
analysis for major arctiine lineages optimized on the ML 
topology, implemented in the program RASP (Yan etal., 2011). 

Ancestral PA acquisition strategy reconstructions with highest marginal 
probabilities are indicated at each node. Blue = Adults associated with 
PAs; Orange = Adults and larvae associated with PAs (e.g., Nyctemera); 
Purple = Larvae associated with PAs (generalists); Green = Larval 



obligate early polyphagous late (e.g., Haploa); Tan = Larvae associated 

with PAs (specialists). 

doi:1 0.1 371/journal.pone.01 01 975.g005 

acquisition has associated evolutionary costs or advantages. 
However, it is known that at least one Haploa species, H. clymene 
(Brown), uses PA-plant derivatives in courtship behaviors that are 
similar to other PA-plant feeding larvae [73]. Nyctemera comprises 
species that are obligate PA feeders as larvae [2], [61], [69] and 
are also facultative PA feeders as adults (Figure 5; [2]). 

Generalist larval feeding on PA plants could be the ancestral 
condition for the Arctiina clade (Figure 5). This habit appears to 
be widespread across the subtribe as it occurs in all generic groups 
except Virhia; the apparent absence of PA feeding in Vii'bia could 
however, be attributed to scant feeding records for the group. The 
fall webworm moth, Hyphantria cunea, is known to feed on over 
400 species of plants [13], [74], although its associations with PA 
plants are not well documented. The salt marsh moth, Estigmene 
acrea, a well-known vegetable pest species [75-76], is pharmaco- 
phagous as larvae [77]. Species of Spilosoma are highly 
polyphagous as larvae with occasional documentation of PA 
pharmacophagy [78]. Other genera such as Apantesis, Arachnis, 
Pyrrharctia are all considered general feeders that likely faculta- 
tively switch between PA and non-PA host plants [3] . In summary, 
larval pharmacophagy evolves at least twice (Figure 5) in this part 
of the tree, obligate larval PA feeding evolves three times (Figure 5) 
and adult pharmacophagy evolves once in Nyctemera (Figure 5). 

Only adult pharmacophagy has been documented for species in 
Pericopina, Euchromiina and Ctenuchina, but verified larval host 
records are rare. Adult pharmacophagy originates three times with 
cases present in all three subtribes (Figure 5). Within the 
Phaegopterina, adult PA feeding appears to be most common in 
the Halysidota group with adult acquisition documented for 
Leucanopsis, Halysidota and Pseudohemihyalea [2], [5-6], [79- 
80]. Within the Euchromina, adult PA feeding is derived and 
restricted to members of genus Euchromia (Figure 5). Adult species 
in this genus have been observed feeding PA plants in the Solomon 
Islands (Plate 10 [A-C] in reference [11] and feeding on 
Heliotropium plants in East Africa [28]. Adult pharmacophagy is 
reconstructed as the ancestral condition for Ctenuchina with 
occurrences documented for Ctenucha [80-81], Antichloris species 
have been observed pollinating and visiting the seeds and dead 
leaves of PA plants [6] while Hyaleucerea species have been 
documented visiting and feeding at PA plant baits [5]. Adult 
pharmacophagy has been documented in other euchromiines and 
ctenuchines [82] but their inclusion in our analysis was not 
possible due the unavailability of DNA quality material. 

The results from this study indicate that facultative pharmaco- 
phagy on PAs as larvae is a derived condition that likely evolved 
from PA plant specialization (e.g., Uletheisa, Tyria). This pattern 
challenges previous hypotheses that PA specialists arose from a PA 
generalist ancestor [8], [10]. Rather, our topology favors a 
specialist to generalist trend, as hypothesized by KrasnoflF and 
Roelofs [83]. Other authors have suggested that facultative PA 
feeding (or more generally, host switching) is in fact a form of 
specialization allowing larvae to gain nutrition from one or more 
plant sources and self-medicate from another [12]. Larval 
pharmacophagous species are also physiologically adapted to 
deactivate a broad range of PAs [12], [77]. For example, in E. 
acrea (and others), the gustatory system is equipped with 
specialized cells that stimulate PA feeding [84—85]. Singer and 
Bernays [3] discuss additional possible benefits to a 'specialized 
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Our results regarding adult pharmacophagy on PA plants 
confirm earlier studies [.'^-6], [10] that inferred multiple origins of 
this behavior within the subfamily. In examined pharmacopha- 
gous species, this habit is accompanied by apparently unique 
sensory structures on the mouthparts of adults [86]. This recent 
study reported a statistically significant increase in the number of 
putative PA chemoreceptors on the proboscides of PA pharma- 
cophagous adults when compared with non-pharmacophagous 
adults or adults of species that obtain PAs only during the larval 
stages [86] . An intriguing possibility is that micromorphology may 
be used to predict adult pharmacophagous behavior. We consider 
this study to be an important initial step towards unraveling the 
intricate evolutionary- interactions between tiger moths and their 
complex associations with pyrrolizidine alkaloids. 

Supporting Information 

Figure SI Phylogenetic hypothesis for the subfamily 
Arctiinae (Noctuoidea, Erebidae) based on Bayesian 
Inference (BI), along with outgroups. Clades representing 
tribes are colored. Support values (posterior probabilities) are 
shown next to the branches. 
(PDF) 

Table SI Complete list of specimens used in the 
phylogenetic analysis including the collection locality, 
sample IDs, and the gene fragments that were success- 
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